Entropy Balance and Dispersive Oscillations in Lattice Boltzmann Models

We conduct an investigation into the dispersive post-shock oscillations in the entropic lattice-Boltzmann method (ELBM). To this end we use a root finding algorithm to implement the ELBM which displays fast cubic convergence and guaranties the proper sign of dissipation. The resulting simulation on the one-dimensional shock tube shows no benefit in terms of regularization from using the ELBM over the standard LBGK method. We also conduct an experiment investigating of the LBGK method using median filtering at a single point per time step. Here we observe that significant regularization can be achieved.Comment: 18 pages, 4 figures; 13/07/2009 Matlab code added to appendi

Spectral singularities and Bragg scattering in complex crystals

Spectral singularities that spoil the completeness of Bloch-Floquet states may occur in non-Hermitian Hamiltonians with complex periodic potentials. Here an equivalence is established between spectral singularities in complex crystals and secularities that arise in Bragg diffraction patterns. Signatures of spectral singularities in a scattering process with wave packets are elucidated for a PT-symmetric complex crystal.Comment: 6 pages, 5 figures, to be published in Phys. Rev.

A unified hyperbolic formulation for viscous fluids and elastoplastic solids

We discuss a unified flow theory which in a single system of hyperbolic partial differential equations (PDEs) can describe the two main branches of continuum mechanics, fluid dynamics, and solid dynamics. The fundamental difference from the classical continuum models, such as the Navier-Stokes for example, is that the finite length scale of the continuum particles is not ignored but kept in the model in order to semi-explicitly describe the essence of any flows, that is the process of continuum particles rearrangements. To allow the continuum particle rearrangements, we admit the deformability of particle which is described by the distortion field. The ability of media to flow is characterized by the strain dissipation time which is a characteristic time necessary for a continuum particle to rearrange with one of its neighboring particles. It is shown that the continuum particle length scale is intimately connected with the dissipation time. The governing equations are represented by a system of first order hyperbolic PDEs with source terms modeling the dissipation due to particle rearrangements. Numerical examples justifying the reliability of the proposed approach are demonstrated.Comment: 6 figure

Quasinormal ringing of acoustic black holes in Laval nozzles: Numerical simulations

Quasinormal ringing of acoustic black holes in Laval nozzles is discussed. The equation for sounds in a transonic flow is written into a Schr\"{o}dinger-type equation with a potential barrier, and the quasinormal frequencies are calculated semianalytically. From the results of numerical simulations, it is shown that the quasinormal modes are actually excited when the transonic flow is formed or slightly perturbed, as well as in the real black hole case. In an actual experiment, however, the purely-outgoing boundary condition will not be satisfied at late times due to the wave reflection at the end of the apparatus, and a late-time ringing will be expressed as a superposition of "boxed" quasinormal modes. It is shown that the late-time ringing damps more slowly than the ordinary quasinormal ringing, while its central frequency is not greatly different from that of the ordinary one. Using this fact, an efficient way for experimentally detecting the quasinormal ringing of an acoustic black hole is discussed.Comment: 9 pages, 8 figures, accepted for publication in Physical Review

A Multi-dimensional Code for Isothermal Magnetohydrodynamic Flows in Astrophysics

We present a multi-dimensional numerical code to solve isothermal magnetohydrodynamic (IMHD) equations for use in modeling astrophysical flows. First, we have built a one-dimensional code which is based on an explicit finite-difference method on an Eulerian grid, called the total variation diminishing (TVD) scheme. Recipes for building the one-dimensional IMHD code, including the normalized right and left eigenvectors of the IMHD Jacobian matrix, are presented. Then, we have extended the one-dimensional code to a multi-dimensional IMHD code through a Strang-type dimensional splitting. In the multi-dimensional code, an explicit cleaning step has been included to eliminate non-zero $\nabla\cdot B$ at every time step. To estimate the proformance of the code, one- and two-dimensional IMHD shock tube tests, and the decay test of a two-dimensional Alfv\'{e}n wave have been done. As an example of astrophysical applications, we have simulated the nonlinear evolution of the two-dimensional Parker instability under a uniform gravity.Comment: Accepted for publication in ApJ, using aaspp4.sty, 22 text pages with 10 figure

Characteristic form of boost-invariant and cylindrically non-symmetric hydrodynamic equations

It is shown that the boost-invariant and cylindrically non-symmetric hydrodynamic equations for baryon-free matter may be reduced to only two coupled differential equations. In the case where the system exhibits the cross-over phase transition, the standard numerical methods may be applied to solve these equations and the proposed scheme allows for a very convenient analysis of the cylindrically non-symmetric hydrodynamic expansion.Comment: 8 pages, 3 figures, 3 sets of figure

Consistent thermodynamic derivative estimates for tabular equations of state

Numerical simulations of compressible fluid flows require an equation of state (EOS) to relate the thermodynamic variables of density, internal energy, temperature, and pressure. A valid EOS must satisfy the thermodynamic conditions of consistency (derivation from a free energy) and stability (positive sound speed squared). When phase transitions are significant, the EOS is complicated and can only be specified in a table. For tabular EOS's such as SESAME from Los Alamos National Laboratory, the consistency and stability conditions take the form of a differential equation relating the derivatives of pressure and energy as functions of temperature and density, along with positivity constraints. Typical software interfaces to such tables based on polynomial or rational interpolants compute derivatives of pressure and energy and may enforce the stability conditions, but do not enforce the consistency condition and its derivatives. We describe a new type of table interface based on a constrained local least squares regression technique. It is applied to several SESAME EOS's showing how the consistency condition can be satisfied to round-off while computing first and second derivatives with demonstrated second-order convergence. An improvement of 14 orders of magnitude over conventional derivatives is demonstrated, although the new method is apparently two orders of magnitude slower, due to the fact that every evaluation requires solving an 11-dimensional nonlinear system.Comment: 29 pages, 9 figures, 16 references, submitted to Phys Rev

The prismatic Sigma 3 (10-10) twin bounday in alpha-Al2O3 investigated by density functional theory and transmission electron microscopy

The microscopic structure of a prismatic $\Sigma 3$ $(10\bar{1}0)$ twin boundary in \aal2o3 is characterized theoretically by ab-initio local-density-functional theory, and experimentally by spatial-resolution electron energy-loss spectroscopy in a scanning transmission electron microscope (STEM), measuring energy-loss near-edge structures (ELNES) of the oxygen $K$-ionization edge. Theoretically, two distinct microscopic variants for this twin interface with low interface energies are derived and analysed. Experimentally, it is demonstrated that the spatial and energetical resolutions of present high-performance STEM instruments are insufficient to discriminate the subtle differences of the two proposed interface variants. It is predicted that for the currently developed next generation of analytical electron microscopes the prismatic twin interface will provide a promising benchmark case to demonstrate the achievement of ELNES with spatial resolution of individual atom columns

Continuous-time link-based kinematic wave model: formulation, solution existence, and well-posedness

We present a continuous-time link-based kinematic wave model (LKWM) for dynamic traffic networks based on the scalar conservation law model. Derivation of the LKWM involves the variational principle for the Hamilton-Jacobi equation and junction models defined via the notions of demand and supply. We show that the proposed LKWM can be formulated as a system of differential algebraic equations (DAEs), which captures shock formation and propagation, as well as queue spillback. The DAE system, as we show in this paper, is the continuous-time counterpart of the link transmission model. In addition, we present a solution existence theory for the continuous-time network model and investigate continuous dependence of the solution on the initial data, a property known as well-posedness. We test the DAE system extensively on several small and large networks and demonstrate its numerical efficiency.Comment: 39 pages, 14 figures, 2 tables, Transportmetrica B: Transport Dynamics 201

Relativistic MHD Simulations of Jets with Toroidal Magnetic Fields

This paper presents an application of the recent relativistic HLLC approximate Riemann solver by Mignone & Bodo to magnetized flows with vanishing normal component of the magnetic field. The numerical scheme is validated in two dimensions by investigating the propagation of axisymmetric jets with toroidal magnetic fields. The selected jet models show that the HLLC solver yields sharper resolution of contact and shear waves and better convergence properties over the traditional HLL approach.Comment: 12 pages, 5 figure